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^ ; Abstract 

■ Since 1978 superconducting coupled cavities have been proposed as sensitive detector of gravi- 

^ ' tational waves. The interaction of the gravitational wave with the cavity walls, and the resulting 

^ . motion, induces the transition of some electromagnetic energy from an initially excited cavity mode 

'nI" \ to an empty one. The energy transfer is maximum when the frequency of the wave is equal to the 

\ frequency difference of the two cavity modes. In this paper the basic principles of the detector are 

' discussed. The interaction of a gravitational wave with the cavity walls is studied in the proper 

' reference frame of the detector, and the coupling between two electromagnetic normal modes in- 

. duced by the wall motion is analyzed in detail. Noise sources are also considered; in particular 

O , the noise coming from the brownian motion of the cavity walls is analyzed. Some ideas for the 

developement of a realistic detector of gravitational waves are discussed; the outline of a possible 
detector design and its expected sensitivity are also shown. 
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INTRODUCTION 



In a series of papers it was studied how the effects due to the interaction between the 
gravitational and the electromagnetic fields could be used to detect gravitational waves 

0]. The proposed detector exploits the energy transfer induced by the gravitational wave 
between two levels of an electromagnetic resonator, whose frequencies uji and uj2 are both 
much larger than the angular frequency VL of the g.w. and satisfy the resonance conditon \uj2— 
cuil = ^ In the scheme suggested by Bernard et al. the two levels are obtained by coupling 
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^ The interaction between the g.w. and the detector is characterized by a transfer of energy and of angular 
momentum. Since the elicity of the g.w., i.e. the angular momentum along the direction of propagation, 
is 2, it can induce a transition between the two levels provided their angular momenta differ by 2; this can 
be achived by putting the two cavities at right angle or by a suitable polarization of the electromagnetic 
field inside the resonator. 
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two identical high frequency cavities^; the angular frequency ui is the frequency of the level 
symmetrical in the fields of the two cavities, and 002 is that of the antisymmetrical one. The 
frequency difference between the symmetric and the antisymmetric level is determined by the 
coupling, and can be adjusted by a careful resonator design. Since the detector sensitivity is 
proportional to the square of the resonator quality factor, superconducting cavities should 
be used for maximum sensitivity. 

The power transfer between the levels of a resonator made up of two pill-box cavities, 
mounted end-to-end and coupled by a small circular aperture in their common endwall, 
was checked in a series of experiments by Melissinos et al., where the perturbation of the 
resonator volume was induced by a piezoelectric crystal P, Q . Recently the experiment was 
repeated by our group with an improved experimental set-up; we obtained a sensitivity to 
fractional deformations of the resonator length as small as 5i/i ~ lO"^*' Hz~^/^ 0. 

In this paper we shall discuss the mechanism of the interaction of a gravitational wave 
with a detector based on two coupled resonant cavities. In previous works this issue was 
discussed using the concept of a dielectric tensor associated with the gravitational wave . 
The interaction was analyzed in the reference frame where the resonator walls were at rest 
even in presence of a gravitational perturbation. We shall analyze the effect in the proper 
reference frame attached to the detector and we shall therefore consider the interaction 
between the wave and the field stored inside the resonator due to the coupling of the g.w. 
with the mechanical structure of the detector [0. 

The paper will be organized as follows: in section |1| the problem of finding the electro- 
magnetic fields in a closed volume with time-varying boundary conditions is studied, and 
an approximate expression of the normal modes in a perturbed resonator is worked out. In 
section |TT| we shall analyze the interaction of a g.w. with the mechanical structure of the 
detector; we shall see that the transfer of energy between a mechanical and an electromag- 
netic oscillation depends both on the electromagnetic field distribution inside the resonator 
and on the resonator geometry and mechanical properties. In section ^ we shall discuss 
some ideas for the developement of a realistic gravitational wave detector based on spherical 
microwave cavities. Afterwords the coupled equations of motion for the fields in a perturbed 
resonator are worked out and solved. In section VIII the issue of the thermal noise of the 
detector's walls is studied; other noise contribution are also considered. Finally, in the last 
section, the expected sensitivity of some detector configuration is shown and discussed. 



II. ELECTROMAGNETIC FIELD IN A RESONATOR WITH PERTURBED 
BOUNDARIES 

To study the mechanism of the energy transfer between the two levels of an electromag- 
netic resonator perturbed by a gravitational wave we shall follow classic electromagnetic 
theory. We shall make use of the fact that any field configuration inside the resonator can 
be expressed as the superposition of the electromagnetic normal modes of the given resonator 

If no sources are present, the electromagnetic field in vacuum is determined by the equa- 



^ Throughout this paper, we shall call resonator the whole detector made up of two coupled cavities; e.g. 
we shall speak about one resonator composed by two coupled spherical cavities. 
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tions: 

V-^ = (1) 
V-H = (2) 

dH 

VAE + fio— = (3) 

VAH-eo^ = (4) 

As can be easily verified, eqs. are automatically satisfied if the fields satisfy the wave 

equations: 

V^^-i^^O (6) 

with c = (/ioco)^^^^. 

Let us assume that the field is contained in a resonator with perfectly conducting walls. 
If we impose boundary conditions on the fields, the solution of the wave equations (H)-® 
will have an infinite discrete set of normal-mode solutions orthogonal to one another and 
complete, in the sense that any arbitrary field in the resonator can be expressed as a sum 
of these normal modes with suitable amplitudes. The amplitudes of each mode can then be 
used to describe the field in the resonator. 

By the familiar procedure of separation of variables we may assume a solution of eqs. 
(|)-(|) of the form: 

oo 

Eif,t) = J2£n{t)E^if^ (7) 



n=0 



and 



where we have defined: 



and 



H{f,t) = Y.nn{t)Hn{^ (8) 

n=0 

£^{t) = ^J E-EndV (9) 



nnit) = ^o H-H^dV (10) 



V 



where the integrals are performed over the resonator volume. 

We require that at the walls the tangential component of E and the normal component 
of H vanish. With this assumption the functions En{r) and Hn{r) satisfy the equations: 

krrEn = V A i?„ (11) 
= V A i„ (12) 

where /c„ = Un/c is the propagation constant associated with the rith mode. 
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It can be proved that the normal modes En and Hn have orthogonahty properties of the 
form 

En-EmdV = 6nm (13) 
Hn ■ HmdV = 6nm (14) 

For a cubical resonator the functions En{f) and Hn{r) are sin(/c„ ■ r) and cos(/c„ ■ r). For 
other geometries they will be other complete sets of functions. The boundary conditions and 
geometry determine the different modes which are distinguished by the index n. In general 
three numbers are needed to specify a mode; n is an abbreviation for this set of numbers. 

Let us now expand the fields E, H in terms of the orthogonal functions En and Hn'-, when 
we substitute these expansions in Maxwell's equations and equate coefficients, so as to get 
the differential equations satisfied by the various coefficients, we find that equations ([|) and 
(I) are automatically satisfied. From equations (|) and @) we find the following equations 
for the expansion coefficients: 

^^-U;n£nit) = -^nnit) (15) 
at 'sin 

+ COnKnit) = (16) 

We have taken into account the dissipation arising from the finite conductivity of the 
walls introducing the electromagnetic quality factor 

_ UJnfJ'O Jy Hn dV _ Gn_ /..-x 

Rs is the material-dependent surface resistance of the walls, and the geometric factor Gn of 
the nth mode is defined as: 

As can readily be seen, equations (|T5|)-(|l6[) for the field expansion coefficients are decou- 
pled: the modes are independent from one another and behave as simple damped harmonic 
oscillators. 



A. Perturbation of boundaries 

Let us suppose that the resonator's boundary is perturbed so that the eigenvalues and 
eigenfunctions of the perturbed resonator differ but little from those of the original one. 
The perturbation method allows to find the eigenvalues and eigenfunctions of the perturbed 
problem from the knowledge of the original ones. It basically consists in expanding in power 
series of a perturbation parameter a the new modes and frequencies 0: 

i;(f) = En{f^ + aen{f^ + 0{a^) 

#;(f) = Hnif) + a hni^ + Oia^) (19) 

k'n = kn + aKn + 0{a^) 
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In the perturbed resonator we shall have: 

oo 

E{f,t) = Y.£nit)EUf^ (20) 

n=0 

and 

oo 

#(f,t) = ^7Y„(t)#;(r) (21) 

n=0 

with: 

£n{t) = ^jE-E'^dV (22) 

and 

n^{t) = ^ f H-nyv (23) 

where the integrals are now performed over the perturbed volume V. 
The perturbed modes will satisfy the following equations: 

k'J'^ = V A (24) 
k'^H'^ = V A (25) 



When the expansions (|T9|) are inserted into the equations (p^) we obtain a series of equa- 
tions determining the various terms of the expansion (|19|). In the first order approximation, 
the perturbed fields may be written: 



AnmE^ir) (26) 



m=l 



-» 1 °° 

0-hn{r) = --CnnHn{r) + ^ BnmHm{r) (27) 
^ m=l 

'^n = ''^knCnn (28) 

where the sums have to be performed for n ^ m. The expansion coefficients have the form 

h h 

A - " " C ('29^ 

"^""^ ^ 1,2 _ 1.2 (^^) 

with 

Cnm = / {Hn ■ Hm — En ■ Em) dV (31) 

J V 

where V = y — is the (algebraic) difference between the perturbed and the original 
volume.^ 

The calculation of the coupling coefficient Cnm depends on how the resonator is deformed 
by an external force. For this reason in the next section we shall briefly review the study of 
the mechanical behaviour of a body under the influence of an external force. 



The derivation of the eigenmodes and eigenvalues of the perturbed resonator has been made assuming a 
static perturbation; in the following we shall assume that it is also correct for a perturbation which has a 
rate of change much slower than the e.m. field characteristic frequency. 
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III. ANALYSIS OF THE MECHANICAL RESPONSE OF THE DETECTOR 



The interaction of a gravitational wave with the mechanical structure of the detector can 
be studied by means of classical, non-relativistic, linear elasticity theory [0, If M(r) 



denotes the displacement of the mass element at point r, relative to the centre of mass of 
the body in its unperturbed state, and f{r,t) is the volume force density which acts on the 
body, the displacement is the solution of the system of partial differential equations: 

P-Q^ - /^^'^ - + V( V ■ ^) = /(^"'' ^) (32) 

with suitable boundary and initial conditions. p{r) is the mass density of the body and A 
and fi are the material's elastic Lame coefficients. In the following we shall adopt null initial 
conditions: 

«(r,0) = ^(f,0) = (33) 

The expansion theorem [0 states that the displacement of a system in response to an 
applied force is equal to the superposition of the normal modes Cai'f^ of the system^: 



uir 



t) = Y.Ur)qait) (34) 



a=l 

— * 

The normal modes ^a(?^ are the eigen-solutions to 

+ (A + ;u) V(V ■ Q = -iulp^a (35) 

with boundary conditions; here a is an index, or set of indices, labelling the mode of fre- 
quency uJa- The modes are normalized so that 

/ Uf^-Uf^p{^dV = MJ^p (36) 

J Vol 

where is the reduced mass of the a mode. For a homogeneous system = M, where 
M is the mass of the system. 

qa(t) is the generalized coordinate of the a mode, obeying the dynamical equation of 
motion: 

Ut) + ^Ja{t) + u;lq^{t) = ^ (37) 

where an empirical damping term, proportional to the system velocity, has been added; fa{t) 
is the generalized force, given by 

fa{t)=f mt)-U^dV (38) 

J Vol 

The solution of equation (0), satisfying the initial conditions (p3D, can be written in 
term of a Green function integral as 



qa{t) = — ^ /* /a(t') sinuoait - t') cxp -^—^ dt' (39) 

MuJn JO Try 



^ In this section we shall label with greek indices the mechanical normal modes of the system, and with 
latin indices the normal modes of the electromagnetic field stored inside the system. 
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with Ta = 2Qa/uJa IS the amphtude decay time of the system. 

In the frequency domain the asymptotic solution of equation (^) is easily found to be: 

/.(a;)/M 

Substituting into eq. (|3^), we find 



It is clear from eq. (|4l| ) , that if we are interested in the displacement in a narrow frequency 
interval uj ± 5a;, only those modes for which a;^ ~ a; (and 7^ 0), will give a significant 
contribution. 



A. Interaction of a g.w. with the mechanical structure of the detector 

An incoming gravitational wave manifests itself as a tidal force density acting on the 
mechanical structure of the detector. Given the expression of the gravitational force and 
the mechanical properties of the detector, the resulting deformation can be calculated, with 
the aid of the mathematical apparatus outlined in the previous section. 

We are mainly interested in the evaluation of the coupling coefficient Cij(t). Let us note 
that, for small displacements, we can write the integral over the perturbed volume as a 
surface integral in the form (see fig. |l]): 

J V 

[ {Hi ■ Hj - Ei ■ Ej) u{t) ■ dS (42) 
s 

where the integral in the r.h.s of eq. (^2]) is now performed over the unperturbed detector 
boundary. It is worth noting that this integral can be expressed as a superposition of the 
mechanical normal modes of the system. Using the expansion theorem (|3^) we can write: 

C,,(t) = fiH^■ Hj - Ei ■ Ej) u{t) ■ dS = 
J s 

^ g„(t) / (H, ■ Hj - E, ■ 4) ■dS=Y^ q^{t)Ct^ (43) 



where we have defined the time-independent form factor C^^- as: 
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C"= iH,-Hj-E,,-E,)i^-dS (44) 



If the external force couples strongly only to one mechanical mode of the detector (say the 
m), we can write the simplified expression 

C.,(t)=g„.(t)C™ (45) 



^ We remind that the superscript a labels the mechanical normal mode, while the subscripts i and j label 
the electromagnetic modes. 
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In summary, to have an effective coupling between the two electromagnetic modes we 
need: 



that the generalized coordinate qm{t) is different from zero. If the system is initially 
at rest, this is true if and only if the generalized force, is itself different from zero, as 
is shown in eqs. ( pT]) and (|39|); 



2. that the spatial integral in eq. is different from zero. In section ^ we shall discuss 
how this depends on the symmetries of the electromagnetic field and of the perturbed 
volume. 



IV. DETECTOR DESIGN 

In order to build a realistic detector a suitable cavity shape has to be chosen. From quite 
general arguments a detector based on two coupled spherical cavites looks very promising 
(see fig. I). 

In order to approach the interesting frequency range for g.w. detection, the mode splitting 
(i.e. the detection frequency) will be U2 — uji ^ 10 kHz. The internal radius of the spherical 
cavity will he r ^ 100 mm, corresponding to a frequency of the TEqh mode ~ 2 GHz. 
The overall system mass and length will he M ^ 5 kg and L ^ 0.8 m. The choice of 
these frequencies for the resonator and mode splitting will be also useful in order to test the 
feasibility of a detector working at ~ 200 MHz and at a detection frequency of ~ 1 KHz. 

A tuning cell, or a superconducting bellow, will be inserted in the coupling tube between 
the two cavities, allowing to tune the coupling strength (i.e. the detection frequency) in a 
narrow range around the design value. 

From the point of view of the electromagnetic design the spherical cell has the highest 
geometrical factor, and so the highest quality factor, for a given surface resistance. For the 
TEqh mode of a sphere the geometric factor G has a value G ~ 850 Q, while for a standard 
elliptical accelerating cavity the TMqio mode has a value of G ~ 250 Q. Looking at the best 
reported values of quality factor of accelerating cavities, which typically are in the range 
10^°-10^^, we can extrapolate that the quality factor of the TEqh mode of a spherical cavity 
can exceed Q ~ 10^^. 

From the mechanical point of view it is well know that a sphere has the highest interaction 
cross-section with a g.w. and that only a few mechanical modes of the sphere do interact 
with a gravitational perturbation (the quadrupolar ones) [|ri|]. The mechanical design is 



highly simplified if the spherical geometry is used since the deformation of the sphere is 
given by the superposition of just one or two normal modes of vibration and thus can be 
easily modeled. In fact the proposed detector acts essentially as a standard g.w. resonant 
bar detector: the gravitational perturbation interacts with the mechanical structure of the 
resonator, deforming it. The e.m. field stored inside the resonator is affected by the time- 
varying boundary conditions and a small quantity of energy is transferred from the initially 
excited e.m. mode to the initially empty one, provided the g.w. frequency equals the 
frequency difference of the two modes. We emphasize that our detector is sensitive to the 
polarization of the incoming gravitational signal: once the e.m. axis has been chosen inside 
the resonator, a g.w with polarization axes in the direction of the field axis will drive the 
energy transfer between the two modes of the cavity with maximum efficiency. 

Finally the spherical cells can be esily deformed in order to remove the unwanted e.m. 
modes degeneracy and to induce the field polarization suitable for g.w. detection. The 
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interaction between the stored e.m. field and the time-varying boundary conditions is not 
trivial and depends both on how the boundary is deformed by the external perturbation 
and on the spatial distribution of the fields inside the resonator, as shown by the expression 
of the coupling coefficient (eq. (0)). It has been calculated that the optimal field 
spatial distribution is with the field axis of the two cavities orthogonal to each other (see 
fig. Different spatial distributions (e.g. with the field axis along the resonators' axis) 
give a smaller effect or no effect at all. A more detailed discussion of the coupling coefficient 
calculation is done in section |VI[ 

V. FIELD EQUATIONS 

We shall derive the equations of motion for the fields in the coupled system from a general 



hamiltonian formalism |14]. The hamiltonian of the electromagnetic field inside a resonator 



with perfectly conducting walls can be written, in terms of the fields amplitudes: 

H = \j^ (eoi -E + fioH-H) dV (46) 

If we substitute eqs. (0)-(P) in eq. (|46|) and use the orthonormality condition, the hamilto- 
nian for the field becomes: 

H=\T.{^1 + '>€) (47) 
If we define a generalized coordinate X„ as: 

X. = ^ (48) 

and its conjugate momentum as: 

Pn = £n (49) 

the hamiltonian can be written as: 

n 

which is identical to the hamiltonian of an infinite set of uncoupled harmonic oscillators. It 
can be easily verified that the field equations of motion can be derived from this hamiltonian 
by: 



dH 
dH 

which, in terms of the fields become: 

dSry 



-Pn — ^n-^n 



dt 
d'hin 
"df 



Xn = Pn (51) 

= —LOn'Hn 

= ^n^n (52) 



which are identical to eqs. (p!5|)-(p!6|) if wall dissipation is neglected 
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A. Field equations for the perturbed system 



In order to find the equations of motions for tfie fields in tfie perturbed resonator we sliall 
make use of tlie results obtained in the previous sections. Let us consider again an external 
time-dependent perturbation, whose rate of change is much less than the rate of change of 
the fields. The normal modes in the perturbed resonator are given by eqs. (|19|)-(^). If we 
are looking for the fields in a frequency region where only two electromagnetic modes give a 
significant contribution, the hamiltonian of the system, as a function of the perturbed modes 
amplitudes (see eqs. (|2^ ) and (p3D, will be: 



H 



\{8l + nl + sl + nt) (53) 



If we now substitute in the above expression the expansions (|TP|) we obtain the hamiltonian 
written in terms of the unperturbed modes amplitudes: 



9. 



m 



(C]^Hl + Cl^Hl + 2 C]^HiH2^ — q-m fm (54) 



where the hamiltonian of a mechanical harmonic oscillator, coupled to an external, time- 
dependent force, has been included. 

From this hamiltonian the equations of motion for the fields can readily be obtained: 

- UJ2S2 = --^^2 (56) 
Ctt ^2 

^ + cu.Hi = -u,q^ {C^.Hi + Ci^7^2) (57) 

-1 + UJ2'H2 = -002qm (CJ^^l + C^2'H2) (58) 

dq-m Pm /rn\ 

+ Muiqm = fm-fr: (60) 



dt 

where the dissipative terms have been added by hand and where the term which describes 
the back-action effect of the fields on the walls is given by: 

fr: = I {cTini + CT2HI) + CuH,n2 (ei) 
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VI. COUPLING COEFFICIENT CALCULATION 



The explicit calculation of the coupling coefficient is not trivial for an arbitrary 
deformation of the resonator volume and in general can be done only by numerical methods. 

First let us note a general property of the coupling coefficients. Let us consider the 
resonant modes of the two coupled cavities. As previously noted we have the symmetric 
and the antisymmetric mode; for the former the electric field E and the magnetic field H 
in the first cavity are equal to the electric and magnetic field in the second cavity, while 
for the latter E and H in the first cavity are equal respectively to — £^ and —H in the 
second one. From this it follows that in the definition of C21, the integrand expression 
{H2 ■ Hi — E2 ■ El) - where we remind that the subscript 1 indicates the symmetric mode 
and 2 the antisymmetric mode - is odd over the whole detector volume. For this reason 
if the volume perturbation, over which we perform the integration, is symmetric between 
the two cavities, the coupling coefficient vanishes, because the contributions to the integral 
coming from the two cavities, cancel each other. Otherwise, if the volume perturbation 
is antisymmetric (when one cavity shrinks, the other expands) the two contributions are 
added with the same sign, and the couphng coefficient is maximum. This general property 
suggests that we must find a geometrical configuration of our detector such that the volume 
deformation due to a g.w. is antisymmetric for the two cavities. This was pointed out 
already in previous works, where the argument was based on the fact that, since the g.w. 
carries an angular momentum equal to 2, the angular momenta of the fields of the two modes 
should differ by 2. This can be achieved by putting the two cavities at right angle or by a 
suitable polarization of the electromagnetic field inside the resonator. 

These concepts were verified by both analytical and numerical calculations. General 
arguments suggested that for an ideal spherical hollow resonator, excited in the fundamental 
quadrupolar mechanical mode and in the TEqh electromagnetic mode, we should have^: 

/ (Hi-H2)im-dS^0A 

s 

I {Hi-Hi)U-dS = (62) 

I (H2-H2)im-dS^0 

s 

More detailed calculations, made on a realistic model of the coupled spheres, including 
the central coupling cell and the e.m. input and output ports, were made by finite element 
methods. These calculations showed that CJ^ = < 10"^ while ^ C^i ^ 0.2. 

VII. CALCULATION OF THE DETECTOR'S SIGNAL 

As already pointed out in the introduction, the proposed detector exploits the energy 
transfer induced by the gravitational wave between two levels of an electromagnetic res- 
onator, whose frequencies cui and CU2 are both much larger than the angular frequency ft of 
the g.w. and satisfy the resonance conditon a;2 — a;i ~ il. This is an example of a frequency 



We remind that for a TE e.m. mode we have vanishing electric field on the resonator's surface; for this 
reason the electric field plays no role in the coupling coefficient calculation. 



21 — <^12 — 
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converter, i.e. a nonlinear device in which energy is transferred from a reference frequency 
to a different frequency by an external pump signal. This can be viewed as a three-bodies 
interaction (given by the field- wall interaction term in the hamiltonian (p^) which corre- 
sponds to annihilation of quanta at Ui and Um and creation at UJ2 (or vice versa). For this 
reason we could argue that, since for a small perturbation 'Hi{t) and £i{t) will approximately 
be sinusoidal functions at frequency cui, while 7^2 (^) and £2{t) will oscillate at frequency to'2, 
in the hamiltonian only the terms varying as UJ2 — cui will give a significant contribution to 
the interaction. The d-c terms will just give an average deformation of the detector's walls, 
determining a static frequency shift of the resonant modes, while the rapidly fluctuating 
terms at uj2 + oji would practically average to zero. 

We shall now calculate the field that is excited by the boundary pertubation in mode 
2, starting from an initial condition with mode 1 strongly excited in the resonator. To 
simplify the analysis of the system of differential equations (^)-(|60|) we will neglect the 
small perturbation, due to the external force, on the initially excited e.m. mode (mode 
1), and will set Tii ~ cos(u;it) and £1 ~ sin(a;it), with constant amplitude A^.'^ 
Furthermore we shall consider the coupling between two TE modes of a resonator: for these 
modes we have vanishing electric field on the resonator surface. Switching to the complex 
notation,^ we obtain: 

fr^ ^\c2l'^{n2Hl) (63) 

Finally, to further simplify our calculations, we shall choose a resonator geometry and e.m. 
field distribution so that = = 0. We shall see in a following section that this choice 
is always possible. 

With this assumptions, and taking Qi ~ Q2 = Qi we can recast the coupled system of 
equations in the following form: 

+ ^^2 + ^2^2 = -ujlq^Clini (64) 

^m+Q^qm + U^^qm = ' 3 ^^^1 (65) 



We apply the following substitutions: 



7^2 (t) = A2{t)exp{iuj2t) 

qmit) = Q{t)exp{int) (66) 
fm{t) = F{t)exp{tnt) 



Eqs. (B1)-(|B^) now become: 



A2 + aA2 + bA2 = cQ 

Q + dQ + eQ = gA2 + F/M (67) 



where a, b, c, d, e and g are constant coefficients defined by 



^ Actually the amplitude of mode 1 is kept constant by an external rf power source. 

® In the following it is understood that the physical fields are the real parts of the complex quantities. 



12 



IjJr, 



Q 



—ulC^Ai 



d = 2in + 



L0„ 
O'r. 



2 o2 I • 



Qn 



2 M 



We can now fourier transform eqs. (|67D and solve them for A2{uj). We find: 

. / X cF(uj)/M 

[0 — 00"^ + t auj) [e — uj'^ + t duo) — gc 

7-^2 (^) is then given by: 

re 

'H2{t) = -exp{iuj2t) / 

J — ( 



cF{u)/M exp{iujt) 



(b — iu"^ + i atu) (e — iu"^ + i duj) — gc 



doj 



(68) 



(69) 



(70) 



(71) 



For a plane g.w travelling along the z axis the force density, in the proper reference frame 
attached to the detector, has the form: 



f{x,t) = --p{x) [{d\x + d\y), {d\x - d\y), (0) 



d 1 / ii^m)x X - {^m)y V) pix) dV - 



Vol 



«2 / i{^m)x y + {^m)y x) p{x) dV 



Vol 



If a* (t) is given by: 



h\a{t) h\l3{t) 
h\l3{t) -h\a{t) 



(72) 



where a* (t), is the adimensional amplitude of the wave, and The 
generalized force, acting on the m mechanical mode, then has the form 



where a(t) and f3{t) are sinusoidal functions of frequency i7/(27r), then a^- 
we define the effective lengths of our detector as 



1 

'M Jvoi 
1 



M Jvoi 



{{^m)x X - {^m)y V) p{x) dV 
((^m)rr V + x) p{x) dV 



we can write 



(73) 



(74) 



-f2^a*, and if 



(75) 
(76) 
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and 



= 5{u - n) (77) 



or, making use of eq. ( |68|) 



2M 



The electric field amplitude the initially empty mode can be readily obtained from eq. 
I); the average energy stored in mode number 2 is given by: U2 = (l/2)|7i2p = (l/2)|^2p- 



VIII. NOISE ISSUES 

A. Mechanical thermal noise 

Thermal noise is one of the fundamental limits in the measurement of small displacements. 
In particular it is one of the dominant noise sources in resonant-mass detectors of g.w. and 
a major reason that such detectors operate at cryogenic temperatures. Since our detector 
exploits the coupling of the g.w. with the mechanical structure of the resonator, we have to 
carefully study the thermal noise contribution to our output signal. 

We start again from eqs. (|64D-(|65|) taking now the external force fm(t) as a stochastic 



force with constant power spectrum 5*//, given by ||T5[: 



Making the substitutions: 



^ff = — 7] — I'yj 



7^2 (t) = A2(t)exp(icj2t) 

qm{t) = Q{t) exp{t {002 - LJi)t) (80) 



we obtain the following equations: 

A2 + aA2 + bA2 = cQ 



fit) 

Q + dQ + eQ = 9^2 + ^ exp(z {u^ - uj2)t) (81) 



where a, 6, c, d, e and g are defined as in eq. (|68|) with the parameter VL replaced by the 
difference UJ2 — oJi. 

The first equation in (^) can be solved for Q(t), and we are left with one equation for 
the variable A2{t): 

d*A2 ,d^A2 „d'^A2 ,„ ,dA2 , , 

[d + 0.)— 7-^ (6 + e + ad)——: (do + ea)—— + [eg - 60)^2 = 



dt^ ^ ' dt^ ' ' dt^ ' ' dt 

fit) 



^ exp(z (^1 - ^2)^) (82) 
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For a linear system we can immediately write the spectral density of the amplitude A2 

as 
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SAA{u^) = \p{u^)\-'{-^ySjf (83) 
with 

p{iu) = -u'^ + i{a + d)uj^ + (e + 6 + da)uj'^ -i{ea + db)uj + eg - eh (84) 

being the fourier transform of the linear system's impulse response function. 
From the above equations we readily find the spectral densities of 7^2 (^) as: 

Shh{u;) = b(^ - ^2)r^ {u;2C^,A,f ^^^^ (85) 

B. Other noise sources 

1. Master oscillator phase noise 

To operate our device we have to feed microwave power into one resonant mode (say 
mode 1), in order to detect the energy transfer between the full and the initially empty 
mode driven by the external perturbation. 

To feed power into our device we shall use a voltage controlled microwave oscillator 
locked on mode 1, at frequency uJi. The master oscillator phase noise is filtered through the 
resonator linewidth; the power spectral density has the following frequency dependence : 

9 ^ -1 - 4/3P,/(^iQ) 

^"^^^^ - (l/Q)^ + (o;M-o;,M2 (''^ 

where P/ is the power input level and [3 is the coupling coefficient of mode 1 to the output 
load. From the above equation we can estimate the microwave power noise spectral density 
at the detection frequency UJ2: 

^'"'^^'^ - il/Qy + Mco,-co,M' ^ l^^^J ^^^^ 

This figure can be improved if the receiver discriminates the parity of the e.m. field at 
frequency UI2, i.e. if it is sensitive only to the power excited in mode number 2, rejecting all 
contributions coming from mode number 1. In this way mode 1 becomes decoupled from 
the output load and /5 = 0. The experimental set-up, based on the use of two magic-tees 
which accomplishes this issue is discussed in detail in 0. Of course the mode discrimination 
cannot be ideal, and some power leaking from mode 1 to the detector's output will be 
present. Nevertheless our previous work has demonstrated that with a careful tuning of the 
detection electronics we can obtain /3 < 10~^^ 0. 



2. Amplifier noise 

The input Johnson noise of the first amplifier in the detection electronics has to be 
added to the previous contributions to establish the overall noise spectral density. It can be 
described by the frequency independent spectral density |I5| : 



Sjj = ksT X 10^1°) = ksTe, (88) 
where N is the noise figure of the amplifier (in dB) and T the operating temperature. 
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IX. DETECTOR SENSITIVITY 



The detector sensitivity is ultimately determined by the overall effect of the various noise 
sources discussed in section |VIII| and, eventually, by several others. Of course, depending 



on the characteristics of the system and on the experimental set-up, different noise surces 
will become dominant. 

We shall characterize the noise in our detector by a frequency dependent spectral density 
Sn{f), with dimension Hz~^, defined as follows if a sinusoidal g.w. with known phase 



known frequency / and unknown r.m.s. amplitude \/2ho, impinges on the detector, and 
if we try to detect the wave by fourier analyzing the detector output with a bandwidth A/, 
then the amplitude signal-to-noise ratio will be: 

^ ^" (89) 



2 



N (5„(/)A/)V 

We shall also define the minimum detectable wave amplitude (at 90% C.L.) for a periodic 
source with known frequency and phase as: 

hrmn{f) = l.7 1^^^ (90) 

with dimension Hz~^/^. In the following calculations the average pattern function value 
F = 2/5 has been taken W7\. 



Let us focus our attention on the system mentioned in section |I^ based on two spherical 
niobium cavities working at cui ~ u;2 ~ 2 GHz with a stored energy in the initially excited 
symmetric mode of f/i ~ 10 J per cell. This is a small-scale system with an effective length 
of 0.1 m and a typical weigth of 5 kg. The lowest quadrupolar mechanical mode is at Um ~ 4 
kHz. In the following we shall consider an equivalent temperature of the detection electronics 
Te, = 30 K. 

A possible design of the detector uses both the mechanical resonance of the structure, and 
the e.m. resonance. This can be accomplished if the detector is designed in order to have 
the mechanical mode frequency equal to the e.m. modes frequency difference Um ~ 1^2 — t^i. 
In this frequency range, with reasonable values for the system parameters, the dominant 
noise source will be the noise coming from the brownian motion of the detector walls. The 
expected sensitivity of the detector for co'2 — t^i = uj^ = 4 kHz is shown in figure ^. In figure 
^ the separate contribution of the noise sources - discussed in section |V111| - to the overall 



noise spectral density is shown. We point out that even if in this case the dominant noise 
source is the walls thermal motion a lower Tgq would increase the detection bandwith, as 
shown in figure ^ 

Since our detector is based on a double resonant system (the mechanical resonator and 
the electromagnetic resonator) it can be operated also for frequencies U2 — uji 7^ Um- At 
frequencies 002 — u^i < 1 kHz the master oscillator phase noise will, in general, be dominant 
(see sec. |V111B 1| ), while at frequencies — uji > 10 kHz the noise coming from the 



detection electronics will dominate (at least for T^g ~ 30 K), as shown in figure |^. The 
expected sensitivity of the detector for 002 — uJi = 10 kHz is shown in figure ^. 

In order to work at frequencies UJ2—UJ1 < 1 kHz a large-scale system has to be developed. 
A possible design could be based on two spherical cavities working at 500 MHz, 

with a;2 — cui ~ 1 kHz. This system could have a stored energy of Ui ~ 800 J per cell, an 
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effective length of 0.4 m and a typical weigth of 300 kg. With a rather optimistic (but not 
unrealistic) choice of system parameters one could obtain the sensitivity shown in figure |^. 
We point out that in this figure an electronics equivalent temperature of T^q ~ 1 K has been 
used; also in this case lowering Tf.q corresponds to an increase of the detection bandwidth 
(see fig. |TT|). 

The large-scale system could also be used at higher frequencies; in this case a good 
sensitivity can be achieved in a narrow detection bandwidth (see fig. P^). 

X. CONCLUSIONS 

A first prototype of the detector has been built and successfully tested 0. A detec- 
tor based on two coupled spherical cavities has been designed and preliminar mechanical 
and electromagnetic tests are being made on normal conducting prototypes. The planned 
timeline is as follows: 

• In 2002 a bulk niobium detector (coupled spherical cavities, uo = 2 GHz, uj2 — uji = 10 
kHz, fixed coupling) will be built at CERN; 

• In 2003 a variable coupling detector will be built and tested. 
In the meantime several open problems must be addressed: 

• The mechanical quality factor of the detector has to be maximized in order to suppress 
the noise coming from the brownian motion of the detector walls. Since mechanical 
dissipations arise from materials intrinsic losses and from the coupling of the system 
to the external environment, materials with low intrinsic losses must be used for the 
construction of the detector and the design of a suitable suspension system has to be 
done carefully. 

• The requirement of an high mechanical quality factor has to be matched with the 
requirement of high electromagnetic quality factor. This can be accomplished by the 
use of bulk niobium, which, at low temperatures, has low intrinsic losses both me- 
chanical and electromagnetic, or by the use of a niobium thin film deposited on a high 
mechanical quality factor substrate. Both tecniques present in principle advantages 
and drawbacks. Several prototypes of single-cell, seamles, copper spherical cavities 
have been built at INFN-LNL by E. Palmieri and will be sputter-coated and tested 
at CERN to check the quality of noibium films deposited on spherical substrates. 

• A cryogenic system with a cooling power of ~ 5 W at T 1.8 K and P = 1 bar has to 
be designed. The contribution of the cryogenic system to the noise has to be studied 
carefully. 

• The readout electronics has to be optimized. The use of a low noise transducer, 
possibly based on the SQUID technology, has to be investigated. 
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If experimental results will be encouraging, by the end of 2003 a proposal for the con- 
struction of a g.w detector, based on superconducting rf cavities could be considered. 
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integration volume ■T^ 




perturbed shell 

FIG. 1: Schematic view of the deformed boundary, u is the local displacement vector; dS is a 
vector pointing in the direction of the outer normal of the original surface. 




FIG. 2: Artistic view of the coupled spherical cavities with the central tuning cell. 
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FIG. 3: Electric field magnitude of the TEqh mode. Note the alignment of the field axis. 
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FIG. 5: Separate contribution of various noise sources to small-scale system sensitivity {cOm ^ 
- wi = 4 kHz, Q = 10^°, Qm = 10^, T = 1.8 K, Teg = 30 K). As can be seen in this case the 
sensitivity is determined by the brownian motion of the walls while the deteciton bandwitdh is 
limited by the amplifier noise. 
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FIG. 7: Separate contribution of various noise sources to small-scale system sensitivity {ujm = 4 
kHz, uj2-0Ji = 10 kHz, Q = 10^°, Qm = 10^ T = 1.8 K, T^q = 30 K). As can readily be seen in 
this case the both the sensitivity and the detection bandwidth are limited by the amplifier noise. 
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FIG. 9: Calculated large-scale system sensitivity for a periodic source {ujm ~ tiJ2 — "^i ~ 1 kHz, 
Q = IQio, Qm = 10^, T = 1.8 K, T^g = 1 K). 
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FIG. 10: Separate contribution of various noise sources to large-scale system sensitivity {um ~ 
a;2 - cui 1 kHz, Q = 10^°, Qm = 10^ T = 1.8 K, Teg = 1 K). Note that here the sensitivity is 
limited by the brownian noise while the detection bandwidth is set by the master oscillator phase 
noise. 
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FIG. 11: Detection bandwidth vs. T^q for large-scale system ~ W2 — "^i ~ 1 kHz, Q = 10^", 
= 106, r = 1.8 K). 
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FIG. 12: Calculated large-scale system sensitivity for a periodic source {um ~ 1 kHz, uj2 — ^i = 10 
kHz, Q = 10^0, = 10^, T = 1.8 K, T^^ = 1 K). 
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